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Abstract 

We first give an (5(2"/'^) quantum algoritlim for tlie 0-1 Knapsack problem witii n variables. 
More generally, for 0-1 Integer Linear Programs with n variables and d inequalities we give an 
0(2"/'^n'') quantum algorithm. For d = o{n/ logn) this running time is bounded by 0(2"'^^/'^+'') 
for every e > and in particular it is better than the 0(2"/^) upper bound for general quantum 
search. 

To investigate whether better algorithms for these NP-hard problems are possible, we for- 
mulate a symmetric claw problem corresponding to 0-1 Knapsack and study its quantum query 
complexity. For the symmetric claw problem we establish a lower bound of 0(2"/^) for its 
quantum query complexity. We have an 0(2"/'^) upper bound given by essentially the same 
quantum algorithm that works for Knapsack. 

Additionally, we consider CNF satisfiability of CNF formulas F with no restrictions on clause 
size, but with the number of clauses in F bounded by cn for a constant c, where n is the number 
of variables. We give a 2(1^")"/^ quantum algorithm for satisfiability in this case, where a is a 
constant depending on c. 



1 Introduction 

The goal of the present paper is to explore the possibility of developing quantum algorithms for 
different NP-hard problems that are faster than what we get by a direct application of Grover's 
search algorithm This can be seen as part of an active research theme in quantum computing: 
discovering tasks for which there are quantum algorithms that are significantly faster than classical 
algorithms (c.f. @, §]). The difficulty in the area is that there are as yet only a few known techniques 
for designing quantum algorithms. The two central methods are essentially from Shor's factoring 
algorithm and Grover's search algorithm 0]. Other known algorithms are essentially based on 
these methods. 

In this paper, we consider quantum algorithms for certain NP-hard problems which have a 
special divide-and-conquer structure that can be exploited to design significantly faster quantum 
algorithms than what a direct application of Grover's search algorithm would yield. For example, 
we give an 0(2"/^) quantum algorithm for the 0-1 Knapsack problem with n variables. This is 
a consequence of a more general result for 0-1 Integer Linear Programs with n variables and d 
inequalities for which we give an 0(2"/^n'^) quantum algorithm. For d = o(n/logn) this running 
time is bounded by 0(2"^^/^+'^)) for every e > and in particular it is significantly better than the 



0(2""/^) time bound if we directly apply Grover's search algorithm. Our algorithms are, of course, 
based on Grover's search algorithm and the powerful method of amplitude amplification 0, |5|. 

The next question we address is whether faster quantum algorithms for the above NP-hard 
problems are possible. To study this, we formulate a new claw-like problem (see [Q] for claw 
problems) which we term as the symmetric claw problem and study its quantum query complexity. 
The symmetric claw problem essentially captures the structure of 0-1 Knapsack and we study its 
complexity in the quantum query model j^, ^. For the symmetric claw problem we are able to 
prove a lower bound of 0(2"/^) in the quantum query model using Ambainis' method [^. The 
problem also has an 0(2"/^) upper bound given by essentially the same quantum algorithm that 
works for 0-1 Knapsack. 

Finally, we consider CNF satisfiability for the case of CNF formulas F with number of clauses 
linearly bounded in the number of variables. Note that we do not make any assumptions on the 
size of clauses in F. More precisely, let CNFc denote CNF formulas F such that the number of 
clauses in F is bounded by cn for the given constant c, where n is the number of variables. We 
give a 2^^"'^)"'/^ quantum algorithm for satisfiability of inputs from CNFc, where a is a constant 
depending on c. 

1.1 Preliminaries 

The set {0, 1}" denotes the set of binary strings of length n, and for a positive integer N we denote 
the set {1,2,... ,A^} by [A^]. 

We design new quantum algorithms for the following NP-hard problems. 

0-1 Knapsack 

Input: A list of positive integers ci, C2, . . . , c„ and a positive integer K. 
Problem: Is there a subset S C [n] such that X^ieS ~ ^■ 

Actually, we consider is 0-1 Integer Linear Programs in general. 
0-1 ILP 

Input: Integers Ojj and bi, 1 < j < n and 1 < i < d. 
Problem: Is the following set of linear inequalities feasible? 

n 

aijXj < bi, G {0, 1}, I < j < n 1 < i < d. 

i=i 

We use standard definitions from quantum computing from In particular, we use the 

quantum query model (c.f. |Q, ^). In this model the input is a function / : [A^] — [M] and 
the values of / is accessed by oracle queries. The complexity of computing some property of / is 
measured as the number of oracle queries. A quantum computation with T queries can be seen as 
a sequence of unitary transformations: 

Uo^O ^Ui^ ■■■O ^Ut, 

where the unitary transform O implements the oracle access to / and the Ui are arbitrary unitary 
transformations which do not depend on the input /. 

A central idea from quantum computing, which is a generalized form of Grover's search algo- 
rithm that we use throughout the paper is amplitude amplification |^: Essentially, if we have a 
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quantum algorithm A with success probabihty p then the success probabihty can be amphfied to 
a constant by 0{^/l/p) calls to A and A~^. 

A range tree (see |^] for details) is a data structure for storing a set S oi N elements 
{xi, X2, • • • , xat}, each of which is a d-tuple of integers xi = (xji, Xj2, • • • , a^id)- The data struc- 
ture can be be built in time N{log'^^^N) and needs space N{log'^^^N). The specific property of 
interest that we require is that we can process range queries on a range tree in time log'^A^ to 
retrieve one element, if it exists, of the set S that satisfies the query bounds. More precisely, 
given a d-tuple of real numbers (01,02, . . . ,0^), in time 0{log'^ N) we can search for an element 
Xi = {xii,Xi2, . . . , Xid) in the range tree such that Xij < Oj, 1 < i < d. 



2 New Quantum Algorithms for 0-1 Integer Linear Programs 

Theorem 1 There is an 0(2"/^n'^) quantum algorithm with constant error probability that solves 
0-1 integer linear programs with n variables and d inequalities. 

Proof. Let the input instance be 

n 

aijXj < bi, 1 < i < d, Xj € {0, 1}, 1 < j < n. 

The goal is to find a feasible solution. We give a stepwise description of the quantum algorithm. 

Step 1. Partition the variables into two sets A = {xi, . . . , x^/^} and B = {x^/^^i, . . . , Xn}, where 
note that |^| = n/3 and \B\ = 2n/3. 

Step 2. For each of the 2"/'^ 0-1 assignments / to the variables in A define the d-tuple yj = 
{yi{l). . . . ,yi{d)) of integers, where 

n/3 

i=i 

where Ij is the value of variable Xj in assignment I. 

Step 3. Let N = 2"'^^. The set X = {yj | / is a 0-1 assignment to A} is of size N. In time 
0{N log'^^^ N) build a range tree (c.f. Theorem 2.11]) to store the set X. The range tree 
has size 0(A^log'^~^ N) and range queries can be processed in time log'^ A^. This completes 
the preprocessing phase. 

Step 4. Using the Hadamard transform prepare the uniform superposition over the set T of 0-1 
assignments to the 2n/3 variables in B: 

a6T 



Step 5. Define a unitary transform U as the standard reversible implementation of the following 
classical subroutine / given u € T as input: 

(i) Compute Zj = Yl]=n/3+i o-ij^j, 1 < i < d. 
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(ii) Let Zi = hi — Zi ioi 1 < i < d, giving a (i-tuple of integers (ii, . . . , i^). 

(iii) Search in the range tree for a yj such that < . . . ,yi{d) < z^. Notice that this 
is a range query and it can be processed in time log'^ N in the range tree date structure 
g, Theorem 2.11]. 

(iv) If such a tuple yj is in X then f{u) = 1. Otherwise, f{u) = 0. 

Step 6. With the initial state as the uniform superposition and using the unitary transform [/, 
apply Grover's search algorithm to search for n E T such that f{u) = 1. 

It follows from the well-known analysis of Grover's algorithm and amplitude amplification |^ 
that the above algorithm has running time 0(2"/^ log*^ A^) = 0(2"/^n''), with constant success 
probability. ■ 

Corollary 2 There is an 0(2"/^) time quantum algorithms with constant success probability for 
0-1 Knapsack. Additionally, any NP-hard optimization problem that takes the form of a 0-1 Integer 
Linear Program with constant number of constraints has an 0(2"/^) time quantum algorithms with 
constant success probability. 

Proof. As an instance of 0-1 Knapsack consists of one equation, it can be expressed as a 0-1 
integer linear program with two inequalities. Hence the algorithm of Theorem |^ yields the claimed 
quantum algorithm. 

Consider any optimization problem that has a linear optimality function and the constraints 
can be expressed as a 0-1 integer linear program with d constraints. Using binary search we can 
reduce the optimization problem to feasibility of a 0-1 integer linear program with d+1 constraints, 
which can be solved using the algorithm of Theorem |^ with running time 0(2"^/^) for constant d. 
Thus, the corresponding optimization problem can also be solved in time 0(2"/^) with constant 
success probability. H 

Remark. Interestingly, if the 0-1 ILP is of the following form: 

minimize/maximize ^iCiXi 

s.t. X]j=i '^ij^j = h, 1 < 2 < d, Xj G {0, 1}, I < j < n, 

then it is easy to see that essentially the same quantum algorithm presented in Theorem |l| solves 
this optimization problem in time 0(2"/^), independent of the number of equations d. The reason 
is that we do not have to maintain a range tree data structure as we do not have to process range 
queries. Just a sorted list would suffice as we only need to make exact queries in this case. 

The above kind of 0-1 ILP is referred to as the 0-1 Group Problem in Nemhauser and Wolsey's 
book Q. An interesting instance of this problem is the following NP-hard problem: given a CNF 
formula F (no restrictions on clause size), search for a satisfying assignment that satisfies exactly 
one variable in each clause. This is clearly an instance of the above 0-1 ILP. Thus we have the 
following corollary. 

Corollary 3 There is an 0(2"/'^) quantum algorithm that takes a CNF formula F with n variables 
as input (no restrictions on clause size or number of clauses) and searches, with constant success 
probability, for a satisfying assignment that satisfies exactly one variable in each clause. 
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3 The Symmetric Claw Problem 



In order to study how far we can exploit this idea of dividing the input, we examine in this section 
a black-box version of the 0-1 Knapsack problem which we term as the Symmetric Claw problem. 

For a pair of strings x,y G {0, 1}" let x\y denote the substring of x obtained by projecting it on 
the positions where y is 1. Likewise, let x\y denote the substring of x obtained by projecting it on 
the positions where y is 0. 

Definition 4 The input to the symmetric claw problem is a function P : {0, 1}" x {0, 1}" 

{0, 1}™" X {0, l}*" where we write P{x,y) = {Pi{x,y), P2{x,y)) for every x,y G {0,1}'^. Additionally, 
the input function P is such that it fulfills the following promise symmetry conditions: 

1. For X € {0, 1}" , if Pi{x, y) = P2{x, y) for some y € {0, 1}" then Pi{x, y) = P2{x, y) for every 
ye {0,1}". 

2. For any X, x' ,y' £ {0,1}^ if x\y = x'\y then Pi{x,y) = Pi{x' ,y). Likewise, if x\y = x'\y then 
P2{x,y) = P2{x',y). 

Given such an input P the symmetric claw problem is to find an x € {0, 1}"", if it exists, such 
that Pi{x, y) = P2{x, y) for some y £ {0, 1}"", by querying P on {x, y) G {0, 1}" x {0, 1}". 

Notice that the 0-1 Knapsack problem is an instance of the symmetric claw problem. 

Remark. Consider a simpler version of the 0-1 Knapsack problem, where we seek a 0-1 assignment 
to variables . . . , x„ such that Yl^=i '^i^i — 0> given integers ai, . . . , o„ as input. 

For each y € {0, 1}" define a partition of the variable set ^ = {i € [n] \ yi = 1} and B = {i £ 
[n] I yi = 0}. Define Pi{x, y) = Y^iaA "'i^i ^^^i^^ v) = - YlieB 

It is easy to check that Pi and P2 fulfill conditions of Definition ^ making this version of 0-1 
Knapsack an instance of the symmetric claw problem. 

We first establish the upper bound result for symmetric claw problem. It is the same algorithm 
we described in the previous section. 

Theorem 5 There is an 0(2"/^) quantum algorithm with constant error probability for the sym- 
metric claw problem in the quantum query model. 

Proof. The algorithm is along exactly the same lines as the algorithm in Theorem |l] for 0-1 ILP. 
In fact it is an easier version of that algorithm as we do not require range trees here. 

We fix y G {0, 1}" to be such that yi = 1 for 1 < i < n/3 and yi = for n/3 -|- 1 < z < n. 

Let Xa = {x e {0, 1}" I Xi = for n/3 + 1 < i < n} and Xb = {x e {0, 1}" \ Xi = for 
1 < i < n/3}. 

Compute and sort the list {Pi{x, y) \ x G Xa} of size 2"/^ using a classical sorting algorithm in 
time 0(2"/^). This is the preprocessing phase. 

Prepare a superposition of x G Xb and search for P2{x,y) in the sorted list using a classical 
binary search implemented by a unitary transform U defined as: For x G Xb, U\x) = 1 if P2{x,y) 
occurs in the sorted list and U\x) = otherwise. 

Now, using Grover's search we can find x G Xb such that Pi{x' , y) = ^2(2;, y) for some x' G Xa. 
At this point, using the properties of Pi and P2 we can put together a complete 0-1 assignment x" 
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from X and x' such that Pi{x",y) = P2{x",y). We can obtain x" by concatenating the first n/3 
bits from x' and the last 2n/3 bits from x. ■ 

We now turn to the question of lower bounds. Using the technique of Ambainis Q we establish 
a lower bound of 17(2"/^) for the quantum query complexity of the symmetric claw problem. We 
first state the result in a form that we need for our setting. 

Theorem 6 H] Let T he the set of input functions f , where f : [M] [M], and let (j) : ^ X x X 
be the function which we wish to compute. Let X,Y be two subsets of T such that (pif) ^ 4'ig) for 
all f e X and g ^Y. Let R X x Y such that: 

(i) For every f X there are at least m different g (zY such that {f,g) G R. 

(a) For every g (zY there are at least m' different f & X such that {f,g) G R. 

(Hi) For every f X and x £ [N], there are at most I different g (z Y such that {f,g) G R such 
that f{x) / g{x). 

(iv) For every g (zY and x £ [N], there are at most I' different f £ X such that {f,g) G R such 
that f{x) + g{x). 

Any quantum algorithm that evaluates <p with constant success probability must make Q,{\J^^-) 
queries to the input function in the quantum query model. 

Now we are ready to prove our lower bound result of this section. 

Theorem 7 Any quantum algorithm that solves the symmetric claw problem with constant success 
probability needs 0(2"/^) queries in the quantum query model. 

Proof. We will prove the claimed lower bound by applying Theorem ^. The inputs consist 
of the class T of functions P : {0,1}" x {0,1}" {0,1}" x {0,1}", where we write P{x,y) = 
{Pi{x,y), P2{x,y)), and every P in ^ fulfills the conditions of the Definition ^. The goal of a 
quantum algorithm (p for the symmetric claw problem is to find x G {0, 1}" such that Pi{x,y) = 
P2{x, y) for some (and therefore, by the promise of Definition ^, every y). Let NO denote the subset 
of "no" instances of J- and YES the subset of "yes" instances. Define X C NO as: 

X = {P G ^ I Pi{x,y) ^ P2{x,y) Vx,y G {0, 1}", Pi(x, y) / Pi{x',y) if x\y / x'\y and 
P2{x,y) / P2{x',y) if x\yj^ x'\y}. 

We will define the subset Y of YES by first defining a relation P C X x YES and then setting 
y = {P' G YES I (P, P') G R for some P G X}. We now define the relation R. 

Given P G X, for every x G {0,1}" we define a distinct P' G YES and include {P,P') in R. 
Given P £ X and x G {0, 1}", P' is obtained from P as follows: 

• If Z]r=i y« — ^/^ then define P[{x,y) = Pi{x,y), and P2(x,y) = Pi{x,y). Furthermore, for 
every x' G {0.1}" such that x'\y = x\y define P2{x',y) = Pi{x,y). For such x' , notice that by 
definition of X we have P{{x',y) ^ Pi{x,y). 

• If J2'i=iyi < then define P2(x,y) = P2{x,y), and P{{x,y) = P2{x,y). Furthermore, for 
every x' G {0.1}" such that x'\y = x\y define P{{x',y) = P2{x,y). For such x', notice that by 
definition of X we have P2{x',y) / P2{x,y). 
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• For all other pairs {u, v) G {0, 1}" x {0, 1}" define P[{u, v) = Pi{u, v) and P2{u, v) = P2{u, v). 

Let P' be defined as above for a given P S X and x G {0,1}". Notice that by construction 
P' is in YES and, furthermore, x is the unique element of {0, 1}" such that P[{x,y) = P2{x,y) for 
y £ {0, 1}". Thus we define 2" distinct elements P' € YES, one for every x € {0, 1}", such that 
(P, P') G R. Therefore, in the terminology of Theorem ^, m = 2". For the remainder of the proof, 
we denote by the element P' £ YES defined as above from an x € {0, 1}" and a given P £ X. 

We next show that for any P £ X and (x, y) G {0, 1}" x {0, 1}" there are at most 2"/^ elements 
P' £Y such that P{x,y) / P'{x,y) and {P,P') G R. In terms of Theorem | this implies I < 2"/^. 
We argue this in two cases: 

Case 1. Suppose X^^Li Ui ^ ^1/2. Then notice from the definition of R that P!, differs from P at {x, y) 
exactly when x\y = z\y. Thus, there are at most 2"/^ such elements P'^. 

Case 2. Suppose 'Y^^=iyi < n/2. Again, notice from definition of R that P'^ differs from P at {x^y) 
exactly when x\y = z\y. Thus, again, there can be at most 2"/^ such elements P'^. 

Finally, notice that /' < m' always holds. 

We can now apply the lower bound result Theorem Since 

mm' m 2" 
~IV~ - J- 2^' 

The lower bound of ^iJ^) given by Theorem ^ yields $7(2"/^). This completes the proof. 



4 Quantum Algorithm for CNF-SAT 

In this section we use ideas from Section |2| to give a simple quantum search algorithm for CNF-SAT 
when the number of clauses is linearly bounded in the number of variables. More precisely, for any 
constant c > we define the following class of CNF formulas: CNFc = {F | F is in conjunctive 
normal form with n variables and P = Ci A C2 A . . . A Cm where m < cn}. 

Theorem 8 For any constant c > there is an 0(2^^~°^"/^) quantum algorithm for satisfiability 
for inputs F G CNF^ where a is a constant such that a < 1/6 and H{a) < ^. Here, H{a) is the 
entropy function defined as —a logg a — {1 — a) logg(l — a) . 

Proof. Let a < 1/6 be a constant to be fixed later, as claimed in the statement. Let the input 
instance from CNFc be F = Ci A C2 A . . . A Cm from CNFc with variables xi, X2, ■ ■ ■ , Xn- Let 
^ = k, where we assume for simplicity of analysis that k is an integer. Partition the variable set 
{xi, X2, ■ ■ ■ , Xn} into k equal size sets ^1,^42,...,^^, where \Ai\ = an, 1 < i < k. Suppose F 
is satisfiable and a* is some fixed satisfying assignment for F. Denote by Ai the set of variables 
{xi, . . . ,Xn} \ Ai for each i. Denote by b* the partial assignment given by a* when restricted to 
variables in A^. Similarly, let c* denote the partial assignment given by a* when restricted to 
variables in Ai. 

Claim. There is an i : 1 < i < k such that the partial assignment c* satisfies at least (1 — a)m 
clauses of F. 
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Suppose the claim is false. Then Ai satisfies at most (1 — a)m clauses of F for each i. Since 
a* is a satisfying assignment, for each i there is a set Si of more than am clauses satisfied by the 
partial assignment b* such that Si H Sj = $ for i ^ j. This is impossible since there are only m 
clauses in F. The claim follows. 

We now give a step-wise description of the algorithm. The algorithm is a loop with index 
i : 1 < i < k, where in the ith iteration it considers the partition of the variable set into and 
Ai. The algorithm succeeds when it considers a partition Ai and Ai such that c* satisfies at least 
(1 — a)m clauses of the input F (such an index exists by the claim). 

Step 1. Partition the variable set into Ai and Ai. 

Step 2. For 1 < / < am build a sorted table T; consisting of the following set of pairs (u, b): u € {0, 1}™ 
is am m-bit string with exactly / I's and 6 is a truth assignment to variables in Ai such that b 
satisfies each clause Cj for which uj = 1, where I < j < m. All other pairs {u, b) are discarded. 
Then, using a classical sorting algorithm and treating u in (u, b) as an m-bit integer key, sort 
the pairs in Ti in increasing order. This entire step can be done in time 0((_^^)2""). 

Step 3. Using the Hadamard transform prepare a uniform superposition of the truth assignments v 
to variables in A^. Let S be the set of all such truth assignments. Notice that l^l = 2^^"")". 
The superposition is 




Step 4. Define a unitary transform U as the standard reversible implementation of the following 
classical subroutine /: 

(i) Compute the vector u G {0, 1}™ where Uj = 1 if and only if Cj is not satisfied by the 
partial assignment v. 

(ii) If u has more than am I's then f{v) returns value and stop. 

(iii) Otherwise, if u has I I's in it, 1 < Z < am, then do a binary search for u in the table T/. 
If {u,b) is found in Ti for some b then f{v) = 1 else f{v) = 0. Notice that if f{v) = 1 
then V and b together give a satisfying assignment for F. 

Step 5. Now, starting with initial state as \tp) and using the unitary transform U to implement /, 
apply Grover search to find v e S such that f{v) = 1. 

Step 6. Output the satisfying assignment {b,v). 

To argue correctness it suffices to notice by our earlier claim that for some i, Ai and Ai is a 
partition such that in Step 3, f{v) = l for at least the partial assignment v = c* obtained from a*. 

Thus, if F is satisfiable the algorithm will find a satisfying assignment with constant success 
probability in time 0((^'„)2"") + 0(2^^^"^"/^), where the first term is the preprocessing time in 
Step 2, and the second term is the time for the Grover search. 

Since m < cn, an easy analysis yields (^) < 2^(")'^ for large n. Thus, if we choose a < 1/6 
such that (1 — a)/2 > H{a)c + a, then the overall time taken is 0(2(-'^~")"'/^). It suffices to choose 
a such that H{a) < 33- ■ 
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5 Discussion 



The question that arises is whether we can close the gap between the upper and lower bounds for 
the symmetric claw problem. Also, for the 0-1 ILP problem with n variables and d inequalities 
we do not have a better upper bound than 0(2"/^)n'^. To study this issue in the quantum query 
model we define the simultaneous claw problem. 

Definition 9 Let {fi,gi), 1 < i < d be functions where fi : [N] X and gi : [N] X, where X 
is some set. The simultaneous claw problem is to find an x £ [N] such that fi{x) = gi{x) for each 
i = 1,2, 

We can also define a symmetric version of the simultaneous claw problem similar to Definition]^, 
but for simplicity we focus on this definition. It is known [H| that for d = 1 there is an 0{N^^^) 



upper bound and the recent new techniques of Aaronson ||l| followed by Shi's sharpened results ||ll| 
imply an ^}{N'^^^) lower bound for the problem. 

We show here that the simultaneous claw problem has a quantum upper bound of 
0(Ar3/4iQgdjY). On the other hand, we are not able to strengthen the lower bound of 17 (A^^/'^) 
which already holds for a single claw pair {f,g). The following theorem is easy to prove. 

Theorem 10 There is a quantum algorithm that takes as input a collection {fi,gi),l < i < d of 
functions where fi : [N] — > X and gi : [A^] X , where X is some set, which makes d{N^/'^log'^ N) 
queries in the quantum query model and outputs a simultaneous claw x € [A^], if it exists, such that 
fi{x) = gi{x) for each i = 1,2,. . . ,d. 

We recall the r-to-1 collision problem: given as input a function / : [A^] — > X with the promise 
that / is either 1-1 or r-to-1, the problem is to find a pair x 7^ y € [A^] such that f{x) = f{y). The 
recent results of Aaronson O followed by Shi |11] imply an Vl{N'^/^) lower bound for the 2-to-l 



collision problem. 

Analogous to the simultaneous claw problem we can define the simultaneous 2-to-l collision 
problem which is given as input functions fi : [N] X , 1 < i < d, with the promise that either all 
the fi are 1-1 or all are 2-to-l, find a pair x 7^ y S [A^] if it exists such that fi{x) = fi{y) for each i. 
The Vt{N'^/^) lower bound of Shi g clearly holds but we have nothing better than that. On the 



other hand, we have a straightforward algorithm that gives an upper bound for the problem. 

Theorem 11 There is a quantum algorithm that takes as input a collection fi,l < i < d of 
functions fi : [A^] — > X where X is some set, makes d{N^/'^log'^ N) queries in the quantum 
query model and outputs a simultaneous collision x ^ y G [N] such that fi{x) = fi{y) for each 
i = 1,2, . . . ,d, if it exists. 
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